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Abstract 

We present in analytic form the three- loop 0{a 2 s ) correction to the H — > gg partial width 
of the standard-model Higgs boson with intermediate mass Mh <§; 2M t . Its knowledge is 
required because the 0(a s ) correction is so sizeable that the theoretical prediction to this 
order is unlikely to be reliable. For Mh = 100 GeV, the resulting QCD correction factor reads 

1 + (215/12)ai 5) (M H )/7r + 150.419 (a ( s ] (M h )/tt) « 1 + 0.66 + 0.21. The new three-loop 
, correction increases the Higgs-boson hadronic width by an amount of order 1%. 
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Q^ ■ The Higgs boson, H, is the missing link of the standard model (SM) of elementary particle 

physics. Its experimental discovery would eventually solve the longstanding puzzle as to whether 
Qh| nature makes use of the Higgs mechanism of spontaneous symmetry breaking to generate the 

particle masses. So far, direct searches at the CERN Large Electron Positron Collider (LEP) 
have only been able to rule out the mass range Mh < 65.6 GeV at the 95% confidence level 
(CL) Q. On the other hand, exploiting the sensitivity to the Higgs boson via quantum loops, 
a global fit to the latest electroweak precision data predicts Mh = 149 + _}^ GeV together with a 
95% CL upper bound at 550 GeV § 

sa' 
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The coupling of the Higgs boson to a pair of gluons, which is mediated at one loop by virtual 
quarks ||, plays a crucial role in Higgs phenomenology. The Yukawa couplings of the Higgs 
boson to the quark lines being proportional to the respective quark masses, the ggH coupling 
of the SM is essentially generated by the top quark alone. The ggH coupling strength becomes 
independent of the top-quark mass Mt in the limit Mh «C 2Mj. In fact, in extensions of the 
SM by new fermion generations, this property may be exploited by using the ggH coupling as a 
device to count the number of high- mass quarks Q . In contrast to the electroweak p parameter 
0], the ggH coupling is also sensitive to quark isodoublets if they are mass-degenerate. At 
this point, we also wish to remind the reader that, by the Landau- Yang theorem Q, spin-one 
particles such as the photon or the Z boson cannot couple to two real gluons, while spin-zero 
particles such as the Higgs boson do. 

The prospects for the Higgs-boson discovery at the CERN Large Hadron Collider (LHC) 
vastly rely on the gluon-fusion subprocess, gg — > H, which will be the very dominant production 
mechanism over the full Mh range allowed Q. The cross section of inclusive Higgs-boson 
production in proton-proton collisions, pp — ► H + X, is significantly increased, by approximately 
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70% under LHC conditions, by including its leading-order (two-loop) QCD corrections ]?], ||], 
which are intimately related to the ggH coupling. Under such circumstances, the theoretical 
prediction for this extremely relevant observable can by no means be considered to be well under 
control, and it is an urgent matter to compute the next-to-leading-order QCD corrections at 
three loops, since there is no reason to expect them to be negligible. Recently, a first step in this 
direction has been taken by considering the resummation of soft-gluon radiation in pp — > H + X 

1- 

An important ingredient in this complex research programme is the 0(a 2 ) three- loop cor- 
rection to the ggH coupling. Typical Feynman diagrams that contribute in this order are those 
obtained by attaching two virtual gluons to the primary top-quark triangle. There are also 
other classes of diagrams, and they all come in large numbers. The ggH coupling also appears 
as a building block in the theoretical description of the crossed process, H —* gg, which con- 
tributes to the hadronic decay width of the Higgs boson. In the low to intermediate mass range, 



Mh ^ 150 GeV, this decay mode has a branching fraction of up to 7% |LC], 11]. Observing that 
a Higgs boson in this mass range almost exclusively decays to bb pairs, this number may be 
quickly understood by taking the ratio of the H — > gg and H — > bb partial widths in the Born 
approximation, which gives (a s Mn / 'ism^) 2 /27 '. 

The Oioig) correction to the H — > gg decay width was originally derived jl^] in the limit 
Mh <C 2M( by constructing a heavy-top-quark effective Lagrangian and subsequently confirmed 
by a diagrammatic calculation || and via a low-energy theorem (LET) @ in Refs. @, g. 
This correction consists of two-loop contributions connected with gg production and one-loop 
contributions due to ggg and gqq final states, where q stands for the first five quark flavours. 
In contrast to the H — > qq decay with subsequent gluon radiation, in the H — ► gqq diagrams 
of interest here, the qq pair is created through the branching of a virtual gluon, so that these 
contributions survive in the limit of vanishing g-quark mass. In fact, if all quark masses, except 
for Mf, are nullified, the hadronic decay width of the Higgs boson is entirely due to H —* gg and 
the associated higher-order processes under consideration here. Depending on the experimental 
setup, the heavier quarks Q = c,b may be detectable with certain efficiencies. The secondary 
Q quarks from H —* gg —* gQQ will typically be much softer than the primary ones from 
H — > QQ — * gQQ, which may serve as a criterion to distinguish between these two production 
mechanisms. Alternatively, one may attempt to subtract the gQQ contributions from the QCD- 
corrected H — > gg decay width pj. For simplicity, following Refs. M, O], we shall not consider 
such a subtraction for the time being. Futhermore, as in Refs. [||, |l2|] , we shall concentrate on 
the limit Mg <C 2Mt, which is most relevant phenomenologically. Although the LEP1 lower 
bound on Mjj jl| then implies that ni = 5 light quark flavours contribute at the renormalization 
scale fi = Mh, we shall keep ni arbitrary. Thus, the Born result reads 

where Gf is Fermi's constant. The 0(a s ) correction may be included by multiplying Eq. (|l|) 
with @, Q 



K =1 + °^ 



IT 



95 7 /11 1 \ 



(2) 



For fj, = Mh = 100 GeV, this amounts to an increase by about 66%. Since such a sizeable 
correction is unlikely to provide a useful approximation, it is indispensable to go to higher 
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orders. 

The purpose of this letter is to take the next step by extending Eq. (J2J) to 0(a 2 ). To this 
end, we need to calculate three-loop three-point, two-loop four-point, and one-loop five-point 
amplitudes. The contributing final states are gggg, ggqq, qqq'q' , ggg, gqq, gg, and qq. Typical 
diagrams are depicted in Fig. [l]. 




Figure 1: Typical diagrams generating 0{a 2 s ) corrections to T(H — > gg). Bold-faced (dashed) 
lines represent the top quark (Higgs boson). 



Our procedure is similar to that of Ref. [12]. We construct an effective Lagrangian, C e g, by 



integrating out the top quark. This Lagrangian is a linear combination of certain dimension-four 
operators acting in QCD with five quark flavours, while all M% dependence is contained in the 
coefficient functions. We then renormalize this Lagrangian and compute with it the H —* gg 
decay width through 0(af ). For brevity, we do not list here all operators that enter our analysis 
in intermediate steps. Instead, we immediately proceed to the final version of £ e ff, 

c cS = -<L x ^G x J?nc x [0[] . (3) 

Here, [O^] is the renormalized counterpart of the bare operator 0[ = G^", where G aiiv 
is the colour field strength, the superscript denotes bare fields, and primed objects refer to 
the five-flavour effective theory. C\ is the corresponding renormalized coefficient function, which 
carries all Mt dependence. Note that C\ and [0[] are not separately renormalization-group (RG) 
invariant through the order considered, while their product is. From Eq. (||) we may derive a 
general expression for the H — > gg decay width, 

T(H _> gg) = ^Zcflm([0[] [0[]), (4) 

where {[0[] [0[]) is the vacuum polarization of the Higgs field induced by the gluon operator at 
q 2 = Mjj, with q being the external four- momentum. 

In order to cope with the enormous complexity of the problem at hand, we make successive 
use of powerful symbolic manipulation programs. Specifically, we generate the contributing 
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diagrams with the package QGRAF |H] and convert the output to a form that can be used as 
input for the packages MINCER [JO]] and MATAD [|TJ]] , which solve massless and massive three- 
loop integrals, respectively. The cancellation of the ultraviolet singularities, the gauge-parameter 
independence, and the RG invariance serve as strong checks for our calculation. 

We adopt two independent methods to calculate G\. One is based on the LET |l3| and 
naturally extends the analysis of Ref. [^] by one order in a s . This leads us to consider the 
top-quark contributions to the gluon and ghost propagators as well as the gluon-ghost coupling 
through O(af) with all external four-momenta put to zero. Specifically, we need to compute 
189, 25, and 228 three-loop diagrams, respectively. The external Higgs line is then attached 
through differentiation with respect to the top-quark mass according to the LET. From the 
resulting three expressions, C\ is then obtained by solving a linear set of equations [12]. The 



second method is the brute-force calculation of the 657 three-loop three-point diagrams which 
contribute to C\. Both methods lead to the same result, which upon renormalization reads 
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where a s is defined in the MS scheme and Mt is the top-quark pole mass. Since C\ appears as 
an overall factor in £ fj, it also enters the calculation of the gg — > H parton-level cross section 
at next-to-leading order |]]. We should mention that Eq. (||) disagrees with the corresponding 
result recently found in Ref. |J], although the numerical difference is relatively small. 

In fact, Eq. (||) can be obtained from known results via the following all-order generalization 
[17| of the LET: 
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Here a^ = a^ (//) and af should be understood as expressed through a^' (fi) via the decou 
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and rrit(n) being the running top quark mass. Futhermorc, 
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stand for the MS beta- function [19| and the quark anomalous dimension |2C], respectively. Eq. 
(|5|) is then recovered from Eq. (|6j) by expressing mt(fJ>) in terms of M t according to [21]. 

We now turn to the second unknown ingredient in Eq. (||), Im([OjJ [Oj]). In fact, it is 
convenient to calculate {[0±] [0±]) first and then to take the absorptive part of it. There is a 
total of 403 three-loop diagrams to be evaluated. After renormalization, the result is 
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where £ is Riemann's zeta function, with values £(2) = 7r 2 /6 and ("(3) ~ 1.202. 

We are now in a position to find the 0(a 2 ) term of the K factor in Eq. (S). To this end, 
we insert Eqs. (||) and (|l0|) with q 2 = Mjj into the master formula (||) and factor out the Born 

result of Eq. ([!]). In order to get a compact expression, we also eliminate as (fx) in favour of 

S n i)i 



Qs (//) p8[ and choose = Mff. We thus obtain 
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where we have substituted n\ = 5 in the last step. If we also use the measured values Mt 
175 GeV and af\M z ) = 0.118, and assume M H = 100 GeV, we have 
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We observe that the new C(q 2 ) term further increases the well-known 0(a s ) enhancement by 
about one third. If we assume that this trend continues to 0(o? s ) and beyond, then Eq. (11) 
may already be regarded as a useful approximation to the full result. Inclusion of the new 0(a 2 ) 
correction leads to an increase of the Higgs-boson hadronic width by an amount of order 1%. 
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Equation ( |TT| ) may be RG-improved by resumming the terms proportional to \n{Mf /Mfj) 
as described in Ref. fl2]| . This leads to 



K 



1 + 14.938 ° ! "' M °» + 2.978 ^> + 104.499 K'< M "> 



7T 



7T 



+ 44,491 +7.818 



7T 



7T 



7T 



7T 



(13) 



For the Mh values of interest here (65.6 GeV < Mh <C 2Af t ), this amounts to an insignificant 
reduction of the absolute value of K, by at most 0.6%, for Mh = 65.6 GeV. In particular, the 



second line of Eq. (12) remains valid within its accuracy. 



Finally, we wish to mention that the K factor of Eq. (11) also applies to the neutral CP-even 
Higgs bosons of two-Higgs-doublet models such as the minimal supersymmetric extension of the 
standard model, as long as their couplings to gluon pairs are dominantly generated via top-quark 
loops. 
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